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Abstract 

We numerically analyze the evolution of a two-dimensional dilatonic black hole, within the 



o 

CGHS model. We focus our attention on the finite-mass corrections to the universal evaporation 
rate which applies at the large-mass limit. Our numerical results confirm a previous theoretical 
\Q prediction for the first-order (oc 1/M) correction. In addition, our results strongly suggest that the 

next-order (oc 1/M 2 ) term vanishes, and provide a rough estimate for the third-order term. 
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I. INTRODUCTION 



In the semiclassical theory of gravity, macroscopic non-spinning black holes (BHs) emit 
a thermal radiation corresponding to the Hawking temperature Tjj = hc 3 /(8irkBGM) [T]. 
This amounts to an outflux rate E which is strictly proportional to 1/M 2 . 

However, this simple and universal result is expected to hold only at the macroscopic 
limit (M — > oo), and one may anticipate a finite-mass correction. The origin of this correc- 
tion may be understood as follows: To derive the quantum outflux one has to analyze the 
(backward) propagation of the field's modes on the BH background, from future null infin- 
ity (FNI) back to past null infinity (PNI). Hawking's original derivation uses the classical 
Schwarzs child geometry as the background metric (over which the field's modes are propa- 
gated). This is a reasonable approximation as long as the BH is very massive (compared to 
the Planck mass Mp{). However, in principle, one should instead use the self-consistent semi- 
classical geometry as the background BH metric. The smaller the BH mass M, the larger is 
the expected deviation of the semiclassical geometry from the classical Schwarzschild solu- 
tion. Correspondingly, one should expect a finite-mass correction to the universal Hawking 
outflux, which grows with decreasing M. 

It may be of interest to evaluate this finite-mass correction to the semiclassical outflux. 
For example, it has been argued j2] (in the framework of two-dimensional gravity) that 
this correction actually reveals the fundamental non-thermal character of the semiclassical 
outflux. Such deviations from thermality could be highly relevant to the attempts to estimate 
the possible amount of correlations between the emitted particles (that is, the amount of 
"information" encoded in Hawking radiation). However, it is difficult to calculate such 
finite-mass corrections in the realistic four-dimensional (4D) context. The reason is that 
there is no known general expression for the renormalized stress-energy tensor T a p in 4D, 
which makes it hard to construct the semiclassical BH geometry (over which the quantum 
field's modes are to be propagated). 

The situation is much simpler in the two-dimensional (2D) context, however. Callan, 
Giddings, Harvey and Strominger (CGHS) [3J introduced a formalism of 2D gravity in which 
the metric is coupled to a dilaton field and to a large number N of identical massless scalar 
fields. In this 2D framework T a p is known explicitly, allowing one to translate semiclassical 
dynamics into a closed system of partial differential equations (PDEs) |3J. Although the 
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exact solution to these PDEs is not known, certain approximate solutions have been derived 
[I]. Also, it is possible to numerically integrate these PDEs and thereby to explore 2D 
semiclassical dynamics [21 [5] . Among other things, such a numerical integration allows one 
to analyze the rate of evaporation for finite BH mass as well. 

There is a remarkable difference between 2D and 4D classical BHs: Whereas in the latter 
the horizon's surface gravity k scales as 1/M, in 2D it is independent of the BH mass. As 
a consequence, the (large-M) Hawking temperature is constant (i.e. independent of M) in 
the 2D framework, and so is the outflux E [3|. This contrasts with the 4D case, wherein 
E oc M~ 2 . Though, for the same reasons explained above, this 2D constant outflux only 
holds in the macroscopic limit, and one should anticipate finite-mass corrections. 

In Ref. |4| an approximate solution to the CGHS field equations was constructed, accurate 
to first order in N/M. Based on this approximate geometry, it is possible to derive [6] the 
leading-order (oc M~ l ) finite- mass correction to the 2D constant outflux. The corrected 
value was found to be E = K[l/A + c\KjM B + 0(M B 2 )} (in appropriate units; see below). 
Here K = N/12, c\ > is a certain known coefficient [6], and Mb denotes the Bondi mass, 
which is essentially the remaining BH mass (see section 3 for more details). Note that the 
outflux increases with time, because Mb steadily decreases. 

One of the main goals of this paper is to numerically explore this finite-mass correction 
to E, and to verify the aforementioned 0(1/Mb) theoretical prediction. 

An independent numerical integration of the CGHS system has been carried out recently 
by Ashtekar, Pretorius and Ramazanoglu (APR) [2J. They also explored numerically the 
dependence of the outflux on the BH mass. However, their investigation was restricted to 
the range of relatively small masses, where the relevant oc 1/Mb expansion parameter is not 
quite ^> 1 (which makes it harder to interpret E in terms of inverse powers of Mb)- In the 
numerical analysis presented here we significantly increase the BH mass, by a factor 2.5. 
This allows us to carry out a more detailed analysis of E in terms of inverse powers of Mb- 

When we numerically obtained the (time-dependent) value of E, we were striked by its 
remarkable similarity to the above-mentioned first-order corrected theoretical prediction. In 
fact, it was not possible to visually distinguish between the numerical and theoretical curves 
(see Fig. 1 below). This came to us as a surprise, because one may naturally expect to have 
0(1/M B ) corrections as well — and no such corrections can be seen in Fig. 1. This could 
hint that the coefficient of the 0(1/M B ) term in the outflux function E(Mb) (which has not 
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yet been derived analytically) may actually vanish. This unexpected observation motivated 
us to analyze E(Mb) in more detail, in order to get a better insight into the 0(1/M|) term 
(and possibly also into the 0(1/ M^) term). 

Our detailed numerical results strongly support the conjecture that the 0(1/M§) term 
indeed vanishes, and also provide a crude estimate of the (non-vanishing) 0(1/M^) term. 



II. THE MODEL AND FIELD EQUATIONS 

The CGHS model [3] consists of a two-dimensional metric g a p coupled to a dilaton and 
to a large number N ^> 1 of identical massless scalar fields We express the metric in the 
double- null form, namely ds 2 = —e 2p dudv. The action then reads 

N 

' f, .. f.- .. 4- , 

12' 



i r r i - 

- / dudv e-^(-2p, uv + 40, u 0,„ - A 2 e 2 - - V f ijU h, v + 



N 

PlU PlV 



(2.1) 



The last term in the action expresses the semiclassical effects, derived from the trace anomaly. 

The model also contains a cosmological constant A 2 . Throughout this paper we set A = 1. 
This choice (along with c = G = 1) fully determines the system of units, making all variables 
dimensionless 1 . 

To simplify the field equations we introduce new variables (following Refs. [HE]): R = 
e~ 2 ^ and S = 2(p — <fi). In these variables the model's evolution equations take the form 

R)uv = c Kpyuv 1 (2.2) 
$iuv R Piuv I R-i 

where K = N/12 [and p = (S — \nR)/2 is to be substituted ]. There are also two constraint 
equations: 

R]WW Ryu) ^JW ~^~T WW 0, (2.3) 

where hereafter w stands for either u or v, and T ww is the (ww) component of the renormal- 
ized stress-energy tensor T Q p. From the trace anomaly one obtains [3] (via energy-momentum 



1 Formally it is equivalent to the change of variable p' = p + ln(X), which does not affect the field equations 
otherwise. Note that hereafter we also set /,; = in the field equations, as we are dealing here with the 
evaporation of the BH rather than its formation. 
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conservation) an explicit expression for T ww : 

f ww = K [p, ww -p 2 , w +z w (w)] , (2.4) 

where z w (w) is a certain boundary function (to be determined from the initial conditions). 
Setting K = 0, one recovers the classical evolution equations 

RjUV c ; S )UV 0, (2.5) 

and the constraint equations R, ww = R, w S, w . This set of equations admits a one-parameter 
family of classical solutions [up to gauge transformations of the general form u — > u'(u), 
v — > v'(v) ], which is the two-dimensional analog of the standard Schwarzschild solution. In 
the so-called Eddington coordinates this solution takes the simple form 

R = M + e v ~ u , S = v-u, (2.6) 

where M is a constant. For each M > the solution describes a static BH with mass M. 
Note that in the classical solution 

p(u, v) = ln(l + Me u - V ) (classical). (2.7) 

This implies asymptotic flatness (p — > 0) at both PNI (u — > — oo) and FNI (v — > oo). The 
special case M = (also known as the linear dilaton solution) describes a flat spacetime, 
p(u,v) = 0. 

Next we consider the collapse of a thin shell of mass M > 0. Following CGHS, we 
assume that the shell propagates along an ingoing null line, which we set to be v = 0. At 



the classical level the solution is (2.6) at v > and flat at v < 0. At the semiclassical 



level the geometry is still flat at v < 0. However at v > the classical geometry (2.7) 



is now replaced by a corresponding solution of the semiclassical field equations (2.2, 2.3). 
The characteristic initial data for the semiclassical evolution equations may conveniently be 
prescribed on the collapsing shell and along PNI: These are exactly the same initial data as 
in the classical collapsing-shell solution [S]. 



This flat solution takes the form R = e v (M + e u ), S — v — u. [A gauge transformation u — > vt 
— ln(M + e~ u ) will then bring it to its more standard form R = e v ~ u ,S = v — u' ] 
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III. OUTFLUX AT FNI 



Throughout the rest of the paper we shall use u and v to denote the Eddington-like 
coordinates for the semiclassical solution at v > 0. They are defined by the requirement 



that the solution takes the asymptotic form (2.6) at both (left) PNI and FNI (and in addition 
the collapsing shell is located at v = 0). 

The shell collapse triggers the onset of Hawking radiation, namely a nonvanishing energy 
outflux T(u) = T uu (u,v — > oo) at FNI (it is the same quantity that we denoted by E 



throughout the Introduction) 3 . It may be obtained from Eq. (|2.4l) by 



T(u) = f uu (u,v -> oo) - f uu (u,v = 0) (3.1) 

[a convenient combination which cancels out the boundary function z u (u) ], because presum- 
ably no outflux crosses the collapsing shell. This expression depends on the actual solution 



p{u,v) at v > through the RHS of Eq. (2.4). 



In the macroscopic limit (M — > oo), one may substitute the classical geometry (2.7) for 
p{u,v), obtaining the leading-order expression 

1 



TV ^~ K 

T{u) « — 



(1 + Me u ] 



TfL (zero-order; global). (3.2) 



(0) 



We shall denote the term in squared brackets by F(u) for brevity, and refer to it as the 
transition function, describing the onset of Hawking radiation. It starts from zero at early 
u, and quickly approaches the asymptotic value 1 at large positive u. Thus, the outflux 
quickly approaches the constant asymptotic value 

K 

T w — = T[^ e (zero-order; asymptotic) (3.3) 

This zero-order calculation properly describes the outflux at the large-M limit. However, 
for a finite-mass BH the expression for T is modified because the semiclassical p(u, v) differs 



from the classical solution (2.7). The leading-order semiclassical correction to the classical 
solution for R and S (and hence for p) has been analyzed in Ref. jjj. Based on this, the 
leading-order correction to the outflux was calculated in [6]. The asymptotic (i.e. late-w) 

3 Note that the outflux T is specifically defined here as the (uu) component of T a p with respect to the 
Eddington coordinate u (i.e. the one for which 2g uv — > — 1 at FNI). 
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result was found to be 
T(u) ~ 

where 



K 

T 



1 + 



K 



2M B (u) 



rrila 

1 (1 



late 
(1) 



(lst-order; asymptotic) 



M B (u) = M - J T(u')du' 



(3.4) 



(3.5) 



is the Bondi mass. (Essentially M B (u) denotes the remaining BH mass as "seen" by a 
distant observer [9].) 



Finally, by combining the asymptotic result (3.4) with the transition function F(u), we 



arrive at the global, first-order corrected, expression for the outflux: 

K 



Tu) 



1 + 



rpgl 

i (0) 



rpgl 



(lst-order; global) 



(3.6) 



2M B {u) 

We shall shortly verify this approximate expression for T(u) by comparing it to a numerical 
simulation. 

The global expressions Tf^ and are obviously more effective than their respective late- 
time counterparts T l ^ e and T^' e , as they properly describe the transient stage of the onset of 
Hawking radiation. It should be pointed out, though, that the simpler, late-time asymptotic 
expressions T l ^ e and T^ e also have their advantage: They serve as "universal curves" [2] 
for T(M B ) (at zeroth and lst-order in K/M B , respectively), onto which all evolutions with 
sufficiently large initial mass should converge, regardless of initial conditions. 



IV. NUMERICAL RESULTS 



We numerically explored the semiclassical 2D spacetime of shell collapse, using a second- 
order finite-difference code. The initial mass was taken to be M = 20K. 4 (This is to be 
compared with the initial value M = 8K used by APR [2], and much smaller values M < K 
used in earlier analyses |5j). Full details of the analysis will be presented elsewhere [TT] . 
The domain of integration covers the range of u wherein the Bondi mass decreases from its 

4 The semiclassical CGHS model admits an exact scaling law in which K changes (that is, the number of 
scalar fields changes), and at the same time the various model's variables are rescaled by certain powers 
of K [21 [4]. In particular, the shell's mass M, the Bondi mass Mb, and the outflux T all scale as K. The 
functional dependence of T/K on Mb/K is thus invariant to this rescaling. This scaling law allows one 
to obtain results for all K values from numerical integrations with a single fiducial K value, e.g. K = 1. 
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original value Mb = 20K up to Mb ~ 15K (at larger u — i.e. smaller Mb — numerical 
errors start to grow exponentially). We calculated T(u) by Eq. (3.1 ), 5 and then constructed 



M B (u) from it, through Eq. (3.5). 



Figure 1 displays the numerically-obtained function T(u), compared with the global first- 
order theoretical prediction T?K(u). Both functions T and T?K are plotted against the Bondi 
mass Mb{u). For reference, the simpler (but less precise) approximate expressions T l ^ e , TfL , 
and T$ e are also shown. Remarkably, the numerical curve is visually indistinguishable from 
T?K(u), even in the zoomed figure lb. This observation came to us as a surprise, because 
it would be just natural to expect corrections to TEJMb) of order 0(K/M B ) 2 - The two 
grey curves in Fig. (lb) represent our original expectation for the typical order of magnitude 
of such a putative second-order correction term. 6 The graph shows no signature of such 
a correction term. In particular Fig. lb indicates that if such a second-order term at all 
exists, it must be <C than its naively-expected order of magnitude. 

This observation led us to suspect that perhaps there actually is no second-order finite- 
mass correction to TfL and T l ^ e . Furthermore, it provoked the intriguing possibility that 
perhaps the first-order corrected expression T^* e is the exact expression for the outflux (in 
the asymptotic late-time limit). 

To address these questions, we explored the residual AT = T — Tts with a much higher 
zoom level. Figure 2 displays AT as a function of Mb 7 ■ Note the tiny vertical scale ~ 10~ 6 . 
At such a small scale, the truncation error becomes a significant issue. The numerical 
simulation used a grid step-size of 0.0025 (in both Eddington coordinates u and v). We 
also carried out simulations with larger step sizes 0.005 and 0.01, and verified second-order 
convergence (in a certain range of Mb, displayed in Fig. 2). We then used Richardson 
extrapolation to correct the truncation error in our finest run (step-size 0.0025). It is the 
corrected residual AT which is shown in Fig. 2. 8 

5 In practice we have calculated T uu (u) along several lines of constant v, at sufficiently large v values which 
mimic FNI. These different v — const lines are indistinguishable in the figures below (see also footnote [8j). 

6 The coefficient of this second-order term has not been analytically calculated so far. The gray curves 
in Fig. la represent a fiducial value, obtained by naively extending the (known) zeroth and first-order 
coefficients to the next order as a geometric progression. 

7 For the inspected mass-range, 17.3 < Mb < 18, T^L has already approached its asymptotic form T l ^ e . 
Therefore we can use the residual AT defined above to study the next-order correction term for T late . 

8 Beside the finite step-size (and round-off) there are two other potential sources of errors in our simulations: 
(i) We start the simulation at a finite u = uq (rather than at PNI, uq — ¥ — oo); (ii) we compute T uu at 
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(a) (b) 
1 1 , 0.2585, 1 1 1 1 — 




FIG. 1: Comparison of our numerical results for the outflux T(u) (solid black curve), as a function 
of Mb(u), to certain approximate analytical expressions, (a): The two pink curves represent the 
(global) zeroth-order and first-order theoretical predictions, T?L (solid) and (dashed). The 
two green curves represent the corresponding late-time asymptotic expressions Tfet 6 (solid) and 
rplate (daghed). (b) : A blow-up on Fig. la. The additional light-gray curves represent one's naive 
expectation for the typical magnitude of a putative second-order correction term (see main text 
and footnote therein). Note that the numerical curve is visually indistinguishable from T?L(u), 
even in the zoom level of Fig. lb. (The solid black curve actually displays numerical data taken 
along four different v = const lines throughout the range 24 < v < 30, which are again visually 
indistinguishable . ) 

Two observations emerge from this figure: First of all, a residual AT certainly exists [that 
is, the first-order corrected expression T 1 ^ £ (Mb) is not the exact (late-time, asymptotic) 
expression for T — an issue which was left open in Fig. lb ]. Nevertheless, Fig. 2 is also 
very suggestive that a second-order correction probably does not exist. The residual AT 
is well matched by a term — 0.05(K/A)(K/Mb) 3 = AT( 3 ), though with a relatively large 

finite v = Vfi na i (rather than at FNI, Vfi na i — > oo). The finite-uo error is corrected by the combination 
of two methods: First, we apply first-order weak-field semiclassical correction to the (otherwise classical) 
initial data that we set at u — u |llj ; Second, we run the code with three different Uq values and 
perform a Richardson extrapolation. The finite- v^„ a ; error is handled by trying several different values 
of 24 < Vfi na i< 32 and verifying that our results for AT are unaffected by further increasing Vfi na i. 
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FIG. 2: The residual AT (solid black curve) as a function of the Bondi mass Mb- The matched 
curve AT(3) (oc M^ 3 ) is shown by dashed red line. For comparison we also display two additional 
curves (dashed gray) proportional to M^ 2 and M^ 4 . (The wiggles in the solid numerical curve 
result from the round-off error.) 



numerical uncertainty in the pre-factor, which we estimate as ~ ±25%. (This large relative 
uncertainty is obviously attributed to the tiny overall magnitude of the residual, which is 
smaller than T by a factor ~ 10~ 5 . 9 ) 



V. SUMMARY 



Our numerical results confirm the previous theoretical prediction of the first-order finite- 
mass correction ( 3.4|3.6 ). They further suggest the absence of a second-order correction 



term, and provide a rough estimate for the third-order term. Our final result is the following 



9 The main source of this uncertainty is the finite Uq (see footnote [8| . Since we only used three Uq values, 
it is hard to assess the effectiveness of the associated Richardson extrapolation. 
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third-order approximate expression for T(Mb)' 

( K \ 3 1 

(5.1) 



K 

rplate 

J (3) - J 



K ( K 

1 + 2Mn +C3 [ 2M 



LB 



along with its global counterpart TLs = F(u) T^ e . Here C3 is a dimensionless coefficient 
which we estimate as ~ —0.4 (with about ~ 25% uncertainty). 

It would be desired to find the analogous finite-mass corrections to the Hawking outflux 
from a 4D semiclassical BH, but this is obviously a much harder task. 
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Note added in Proof: After this work was submitted, Ramazanoglu notified us that he 
managed to convert the results of Ref. [2] for the mass-dependent outflux, from the modified 
definitions of outflux and Bondi mass (introduced in Ref. [10J ) to the traditional ones [9]. He 
then carried out an expansion of the (traditional) outflux in inverse powers of (traditional) 



Bondi mass. His results, truncated after third order, agree with our Eq. (5.1) with c 3 
-(5/12) « -0.42 
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